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ABSTRACT 

Measurement of the 21cm hyperfine transition of neutral hydrogen provides a unique 
probe of the epoch of reionization and the Dark Ages. Three major mechanisms are be- 
lieved to dominate the radiation process: emission from neutral hydrogen surrounding 
the ionized bubbles of first galaxies and/or quasars, emission from neutral hydrogen 
inside minihalos, and absorption of diffuse neutral hydrogen against the cosmic mi- 
crowave background. In the present work, by simply combining the existing analytic 
models for the three mechanisms, we investigate the contribution of cross-correlation 
between these three components to the total 21cm angular power spectrum, in the 
sense that neutral hydrogen associated with different radiation processes traces the 
large-scale structures of underlying density perturbations. While the overall 21cm 
power spectrum remains almost unchanged with the inclusion of the cross-correlations, 
the cross-correlation may play a key role in the determination of the 21cm power 
spectrum during the transition of 21cm radiation from emission-dominated phase to 
absorption-dominated phase at redshift z w 20. A significant suppression in the 21cm 
angular power spectrum during this transition is anticipated as the result of negative 
contribution of the cross-correlation between the absorption of diffuse neutral hydro- 
gen and the emission components. Therefore, an accurate prediction of the cosmic 
21cm power spectrum should take the cross-correlation into account especially at the 
transition phase. 
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1 INTRODUCTION 

The 21cm hyperfine transition of neutral hydrogen in the form of either emission or absorption against the cosmic microwave 
background (CMB) provides a unique probe of the history of our universe between the surface of last scattering (z ~ 1000) and 
the end of reionization (z « 6). While the expected brightness temperature of these 21cm signals is two orders of magnitude 
smaller than the CMB temperature and can be easily swamped by extremely strong foreground such as the Milky Way, many 
ambitious radio facilities (e.g. 21CMA0, LOFAE0, MWA@, PAPER0, SKA@, etc) have been constructed and are being planned 
to search for the weak signals primarily from the epoch of reionization (EOR) at 6 < z < 20. Indeed, the Dark Ages and the 
EOR may constitute the last frontier of observational cosmology. 

There have been many efforts - both theoretical and numerical - ov er the past few years aimed at und erstanding various 
phys ical processes in the intergalactic medium (IGM) during the EOR [see Furlanetto. Oh fc BriggslBood and Morales fc Wvithel 
2009 for a review]. In order to make the 21cm signals visible, the spin temperature of neutral hydrogen must differ from the 



1 See |http:/ /cosmo.bao.ac.cn| 

2 See http://www.lofar.orgl 

3 See http://www.haystack.mit.edu/arrays/MWA 

4 See http://astro.berkeley.edu/~dbacker/eor 

5 See http://www.skatclescope.org 
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temperature of the CMB, Tomb, and an effective mechanism exists to couple the spin temperature to the gas temperature. 
Apart from some exotic scenarios such as dark matter particle decay, three mechanisms are believed to play dominant roles in 
the generation of observable 21cm signals fro m the Dark Ages to the EOR: scattering of UV photons of first stars primarily 
t hrou gh the so- called Wouthuysen-F i eld effect JWouthuvsenll952l : lFieldlll959h , heating of IGM by collision mainly within mini- 
halos jlVladau. Meiksin fc Reel 19971 : llliey et alJl2002h. and absorption of uncollapsed IG M against the CM B radiation prior to 
the formation of first cosmic structures ( Scott fc Reeslll990l : lLoeb fe ZaldarriagalEool : Madau et al. 1997; iLewis fc Challinorl 
2007). Among these three mechanisms, the foreground absorption of neutral hydrogen against the CMB dominates the 21cm 
background through the epoch of the Dark Ages down to redshift z « 20, while for 6 < z < 20 the energy release of first 
luminous objects (stars and/or quasars) plays a leading role in generating the 21cm emission in excess of the CMB (e.g. 
Scott fc Reelll990l: iBarkana fc Loebll2004r h During the transition phase or the 'Grey Ages' between the end of the Dark Ages 
and the dawn of reionization at z ~ 20, a combined mechanism may come into effect in the redshifted 21cm background. 

If we decompose the surface brightness temperature of 21cm background from the Dark Ages through the EOR into three 
components: emission from neutral hydrogen surrounding the ionized bubbles of first galaxies and/or quasars T em , emission 
inside minihalos T m h, and absorption of IGM against the CMB T a b, the 21cm power spectrum can be formally determined 
by three auto-correlation terms (T om T em ) + (T m hT m h) + (TkbTkb), an d three cross-correlation terms 2(T em T m h) +2(T cm T a b) + 
2(T m hT a b). In particular, two of the cross-correlations relevant to the absorption have negative amplitudes as a result of 
Ts < Tqmb, where Ts is the spin temperature of neutral hydrogen. Alth ough quite a lot of work has been done in recent litera- 



ture on the theoretical estimate of the three auto-correlation terms 
I 
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2008;: [Lewis fc Challinorll2007l : iMao et al.ll2008t ). it has remained unclear so far to what extent the three cross-correlations con- 
tribute to the total 21cm power spectrum, especially in the Grey Ages. The essentials of three cross-correlations arise from the 
fact that both collapsed and uncollapsed neutral hydrogen may trace underlying gravitational poten tials of large-scale density 
pertur bations, through which their 21cm emission/absorption signals are correlated at large-scales. Furlanetto. Oh fc Briggsl 
1 2006) have actually included the cross-correlation between T cm and T m h in the study of the power spectrum of 21cm fluc- 
tuations, although their work focused mainly on the observational distinction between signatures of minihalos and those of 
IGM, in which massive halos (galaxies) are ignored. In the present study we will concentrate on the evaluation of each of the 
cross-correlations in terms of the 21cm angular power spectrum for different frequencies or redshifts. 

Note, however, that two of the cross-correlations, (T em T a b} and (T om T m h), should be safely treated as a mathematical 
elegance rather than a physically realistic situation through most of the period of the Dark Ages to the EOR. The latter 
requires the co-existence of 21cm emission from 'hot' IGM around first galaxies and 21cm absorption from 'cold' IGM, 
represented by (TemTab), and the co-existence of 21cm emission from 'hot' IGM around first galaxies and 21cm emission from 
'hot' IGM i nside minihalos, represented by (T em T m h), within the same regions. Unless an inhomog eneous heating scenario is 
considered ( Barkana fc Loeb1l2004 ; Pritchard fc FurlanettcfeoO? ; Santos et al.ll2008l : Mesingerll2009h . the two correlation terms 
(T em T a b) and (T em T m h) may become physically meaningless except for the very rapid transition phase or the Grey Ages. In 
this regard, the last cross-correlation between the 21cm emission from minihalos and the 21cm absorption from the diffuse 
IGM, (T m hT a b), is probably the only physically plausible term that is worth investigating during and even before the Grey 
Ages. 

The main difficulty of evaluating accurately the cross-correlations at present within the framework of analytic description 
of the reionization process is perhaps the absence of a unification scheme in which the three mechanisms of generating cosmic 
21cm signals can be incorporated and treated at any redshifts through the Dark Ages to the EOR. The emission scenario 
associated with first galaxies which has been extensively discussed in recent literature applies only to lower redshifts at 
6 < z < 20, while the absorption model of uncollapsed IGM is restricted to higher redshifts beyond z « 20. During the Grey 
Ages when the effect of the cross-correlations between different mechanisms is expected to be visible in the cosmic 21cm 
power spectrum, our theoretical prediction relies on the extrapolation of the current emission and absorption models unless 
we make an attempt at exploring the unification scheme in which the distribution and evolution of both neutral hydrogen and 
its reionization can be correctly described through the entire period of the Dark Ages a nd the EOR including the transition 



phase. It is likely that numerical simulation (e.g. Zahn et al. 2007 : Trac fc Cen 2007 : Iliev et al. 20081 : Santos et al. 2008) 



should be eventually used so that various physical processes can be incorporated in the estimate of the cross-correlations in 
the cosmic 21cm power spectrum. 

Having realized these difficulties and limitations, we still intend to take the three cross-correlations into account in the 
calculation of the cosmic 21cm power spectrum. By simple combining the existing analytic scenarios of generating the cosmic 
21cm signals at the Dark Ages and the EOR, we address the question of what and how large observational effect one may 
expect to see if the three correlation terms are included. Our goal in the present paper is to call the attention of the 21cm 
cosmology community to the possible role of the cross-correlations in the 21cm power spectrum especially at the Grey Ages 
rather than the detailed modelling of the evolution of cosmic neutral hydrogen and reionization process. It is hoped that 
future numerical simulation can provide a more sophisticated treatment of the problem. 
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2 GENERAL FRAMEWORK 



We concentrate our attention on the following three mechanisms only: 21cm emission from the neutral hydrogen surrounding 
the Stromgren spheres of the first galaxies/quasars, 21cm emission from minihalos, and 21cm absorption against the CMB 
from the uncollapsed IGM. If we use the surface brightness temperature to denote the cosmic 21cm radiation, the total 
background can be written as 



T — T cm + T m h + T & \ 



(1) 



We usually work with the two-point correlation function or power spectrum in the Fourier space to extract the considerably 
faint 21cm signals in a statistical manner. This yields 



(TT) — (T cm T om ) + (T m i,T m h) + (T a bT a b) + 2(T em T m h) + 2(T cm T a b) + 2(T , m hT a b). 
or in terms of power spectrum 

P(k) = P cm - cm (fc) + P mh - mh (fc) + Pab-ab(fc) + 2P em _ mh (fc) + 2P em _ ab (fe) + 2P mh _ ab (fc). 



(2) 



(3) 



In what follows we utilize the angular power spectrum C(£) by specifying the observing frequency v (or redshift) instead 
of the 3D power spectrum P(k). Moreover, we adopt a small-angle approximation, which allows us to convert the 3D power 
spectrum P(k) of a perturbation field into the corresponding angular power spectrum C(£) in terms of the Limber's equation. 
If the frequency response function of 21cm experiment is W^ (v) or W ro (r) centered on vo or ro, where r is the comoving 
distance, then two-point angular power spectrum can be obtained by 



W 2 (r) 
dr-^P-P{l/r) 



(4) 



In fact, in the calculation of both auto and cross angular power spectra of the 21cm background, many of the perturbations 
can be related to the underlying matter density fluctuations. Consequently, in most cases P(k) simply represents the power 
spectrum of dark matter, and other coefficients can be absorbed into the window function. For the window function, we will 
take a Gaussian form of 
1 



W (r) - --(r-r )72A/ 

in which the resolution Ar in r is related to the observing frequency bandwidth through 



a — »(s&£)( 
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Mpc. 
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(6) 



We ev aluate the power spectrum of density perturbatio n following the analytical halo approach to include non-linear 
structures 1 Seliak 2000l ; Coorav. Hu fc Miralda-Escude 2000h . The non-linear dark matter power spectrum Pu{k, z) consists 
of a single halo term P^(k, z) plus a clustering term P^(k, z): 
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where ph(M,z,k) is the Fourier transform of dark matter density profile p^(M,z) of a halo M at z, and ~p(z) is the 
mean mass density of the ba ckground universe at z. We adopt the functional form suggested by numerical simulations 
jNavarro. Frenk fc White 1997 ) for Ph{M, z) and fix the concentration parameter in terms of the e mpirical fitting form ula of 
Bullock et al.l (|2002f t . b(M,z s ) is the bias parameter, for which we use the analytic prescription of lMo fc White! (1996). The 
mass function of dark halos d 2 n(M, z)/dMdV is assumed to follow the standard Press-Schechter function 



d 2 n(M,z) 



S c da 



2a 2 



(9) 



dMdV \l i: M a 2 dM 

in which a is the linear theory variance of the mass density fluctuation in sphere of mass M: 

poo 

a 2 = (1/2tt 2 ) / k 2 P^(k)\W R (kR)\ 2 dk, (10) 
Jo 

and Wn(kR) is the Fourier representation of the window function. For a top-hat window function, we have WR.(kR) = 



3[sin(fcP) - kRcos(kR)]/(kR) J . We use the present-day linear matter power spectrum Py[ l (k) given bv lBardeen et al 



(198C 



and a power-law form (oc k" s ) for the primordial matter density fluctuation. Moreover, we adopt the concordance cosmological 



model (ACDM) with the following choice of cosmic parameters: S7m 
and as = 0.84. 



0.27, n A = 0.73, fW = 0.0224, h = 0.71, n s = 1.0 
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Since our attention is concentrated on the evaluation of the 21cm angular power spectrum of each of the cross-correlations 
rather than the detailed modelling of the underlying physical process of 21cm transition, we adopt the existing analytic 
algorithms developed in the recent literature to calculate the cosmic 21cm emission/absorption signal. In what follows we 
summarize very briefly each of the analytic or semi-analytic approaches in the computation of auto-correlation angular power 
spectrum and extend the work to the evaluation of the corresponding cross-correlations. 



3 POWER SPECTRA 

3.1 Auto-correlation: C cm -cm(Q 

We follow essentially the prescription of Zaldarriaga et al.(2004) to estimate the 21cm emission angular power spectrum from 
neutral hydrogen surrounding the ionized bubbles of first galaxies. In their treatment, the reionization process is associated 
with the growth of HII regions around individual galaxies, manifested by a sphere of radius of R. We use an analytic form of 
R — {3 + 0.424[1 — exp[(xn — 0.5)/0.375] 2 ]} /i _1 Mpc to account for the cosmic evolution of R, which roughly matches the 
values adopted by Zaldarriaga et al. (2004) at different redshifts. The overall cosmic evolution of reionization is specified by 
the mean neutral hydrogen fraction, for which a simple analytic expression is assumed: xh(z) = 1/1 + exp[— (z — zo)/Az], 
with a choice of zq — 10 and Az = 0.5 in our calculation. The brightness temperature in excess of Tomb is 



T cm = T (l + S)x H (11) 



To = 23mK 



lib n 
0.02 



in which we have assumed that the spin temperature of neutral hydrogen is much higher than Tcmb , and the peculiar motion 
of the IGM is neglected. The perturbation in T em arises from a combination of fluctuations in both the matter density (1 + <5) 
and the neutral fraction ih, giving rise to the power spectrum 

P(i+s)x H (k) = x H P M (fe) + (x H - x' 2 1 )P XH s(k) + (2h - x%)P Xii {k). (13) 

Namely, the power spectrum of (1 + 8)xu comprises two auto-correlations and one cross-correlation between the two fields. 
An analytic model was used to construct P XIi (k) and P XH s{k) in Zaldarriaga et al. (2004) and a refined model was proposed 
in their subsequent work (Furlanetto et al. 2004a). The observed 21cm brightness temperature in the direction at frequency 
v is thus 

T e <£ s (0,i/) = J drWu(r)T em (9,v). (14) 

Applying the Limber's equation, we can obtain the corresponding angular power spectrum C cm _ m(^) by replacing P(£/r) 
with T 2 P (1+s)xn (£/r) in Eq.(4). 



3.2 Auto-correlation: C„ 



It is argued that collisional excitation of the neutral hydrogen inside minihalos may constitute a non-negligible fraction of 
the cosmic 21cm background, especially at the Grey Ages. Many investigations, based on both semi-analytical treatment 
and numerical simulat ions, have been made towards the modeling of t he radiation process and the signatures of the resultant 



21cm background ( e.g. ICarilli. Gnedin fc Qwenlbooj llliev et al-lbool 120031 : ICiardi et all2006l : iFurlanetto. Oh fc Briggsll2006l : 



Shapiro et al.l l2006') Minihalos are restricted by two mass thresholds: the Jeans mass that sets the low-mass limit on collapsed 



halos, and the upper limit beyo nd which the atomic cooling becomes efficient. In the ACDM cosmological model, the two 
mass limits can be estimated by Shapiro et al. ( 20061 ) 
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(15) 



(16) 



It is unlikely that present and even future radio telescopes/arrays will be able to resolve minihalos because of their 
small size and low surface brightness (< /iJy). We therefore deal with the total 21cm flux of each minihalo, and neglect the 
contribution of the Poisson term in the calculation of power spectrum. In order to simplify the procedure, also guided by the 
result of numerical simulations, we assume that the 21cm flux F of a minihalo is pro portional to its total mass M. We then 
fit the simulated line-integrated flux F versus minihalo mass M of llliev et all |2002h to the following linear relation in the 
Log F-Log M space: 
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logF = logM + E(z) 



(17) 



where the only free parameter E(z), which characteriz es the cosmic evolution of minihalos, is determined a nd extrapol a ted at 
each redshift in terms of Figure 1 of llliev et al.l (|2002l ). Note that the cosmological parameters adopted in llliev et al l (|2002h 
are slightly different from ours. Here we have neglected the possible modification to the above relationship. The overall surface 
brightness temperature of all minihalos is given by 

pMmax 
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(18) 



where ks is the Boltzmann constant, and dV/dzdQ. is the comoving volume per unit redshift per unit solid angle. Note that 
we have converted the total flux into the surface brightness temperature in terms of the Rayleigh- Jeans approximation. The 
angular power spectrum of T m h can be obtained using the Limber's equation 

2 
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2k B J 
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W:(r)f(r)P u (e/r)r 2 dr, 



where fir) = /l I f max F\ 

■' y ' JMmm L 



H 2 (z) 

]b(M,z)dM, and H(z) is the Hubble constant at z. 



(19) 



d^n(M,z) 
dMdV 



3.3 Auto-correlation: C a 



It was formulated long time ago that before the formation of first stars at z « 20, neutral hydrogen may absorb the CMB 
flux since the baryon temperat ure dropped faster than the CMB temperature in the framework of adiabatic cooling of baryon 
with evolution of the universe 1 Field 19591 ). The fluctuation of such a 21cm absorption field arises purely from the underlying 
matter density perturbation, T a b = (dTkb/ 'dS)5. The 21cm absorption surface brightness temperature along the line-of-sight is 



drW v (r)^5(0,v) 
do 



The corresponding angular power spectrum can be derived using the Limber's equation 
= f WS(r) /<ffkb\ 2 



C, 



ab — ab 



P M (£/r)dr. 



(20) 



(21) 



\ dS J 

Finally we can use essentially the approach of Loeb fc Zaldarriaea 1 2004h to calculate the brightness temperature variation 
with respect to S, dT^/dS. 



3.4 Cross-correlation: C cm _ m h(^) 

Neutral hydrogen both surrounding the ionized bubbles of first galaxies/quasars and inside minihalos contributes to the cosmic 
21cm emission. If galaxies/quasars (massive halos) and minihalos trace common gravitational potentials of underlying density 
perturbations at large scales, their 21cm emission signals should be strongly correlated, and the combined effect raises the 
amplitude of the total 21cm power spectrum. Yet, the above scenario is too simplistic because it is implicitly assumed that the 
minihalos may survive within HII regions for relatively long time. This may happen during the Grey Ages when the population 
of first luminous galaxies are rare and evaporation process is still inefficient. Another possibility is that minihalos were not yet 
swallowed by the expanding bubbles of first galaxies/quasars and both of them shared the same density perturbation field. 
The cross-correlation C om -mh(^) eventually vanishes with the evaporation process of minihalos. A combination of Eqs.(ll), 
(12) and (18) yields 



(T em T mh ) = Toil + 5)x H 
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Replacing d Jj^y^ y dMdv ft z)S] and removing the global component, we can easily derive the 3D power spectrum 

\3 



( fe ) = T ° ( 4: ) (1 tfV' [(f + g)P XH s(k) + x s fP tA (k)], (23) 



2k B J H{z) 

in which g = f^™^ F[ d dMdv^ dM. Furthermore, the corresponding angular power spectrum under the small-angle approxi- 
mation is 

C em -^)=To(^y J (1 + Z jf™" {r) [(/ + 9)P» B t{t/r) + x H fP M Wr)]dr. (24) 

Note again, however, that the above expressions apply only to the regimes outsides the ionized bubbles surrounding the first 
galaxies. We use the characteristic radius R in Sec. 3.1 to denote this restriction, i.e. the cross-correlation turns to be zero on 
scales smaller than R. In the calculation of the angular power spectrum, the power vanishes if £ > 2nr/R. 
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3.5 Cross-correlation: C c 



Before the EOR, most of IGM is still neutral, and the mean global fraction of ionized hydrogen at z w 20 is well below 
unity. So, one may have an intuition that there should exist a strong correlation between the 21cm emission generated by 
neutral hydrogen surrounding the ionized bubbles of first galaxies/quasars and the 21cm absorption from the uncollapsed, 
diffuse neutral hydrogen. This arises because both the collapsed and uncollapsed baryons trace the common large-scale density 
perturbations of the universe. Yet, the co-existence of 'hot' IGM around first galaxies and 'cold' IGM within the same region 
vanishes rapidly with the fo rmation of more first generation stars, unless the inhomogeneous nature of Lyq and/o r X-ray 
backgrounds are introduced ( Barkana fc Loeb|[20o3 : IPritchard fc Furlanettdl2007l ; ISantos et alJliooj ; lMesingerll2009l '). Recall 
that first galaxies formed in overdense regions, in which the heating of ambient diffuse IGM by energy release of the first stars 
proceeds more efficiently than that in underdense regions. In other words, the mean neutral fraction of the universe should not 
be used as a unique indicator of the significance of the cross correlation (T em T a b). Indeed, the signature of cross-correlation 
(T em T a b} on the total 21cm cosmic angular power spectrum could only be detected during the transition phase. This arises 
because the population of the first stars even in the overdense regions are rare during this rapid transition process and most 
of the radiation energy is not yet delivered to the uncollapsed IGM in the same regions. Such a cross-correlation makes a 
negative contribution to, and may therefore suppress the total 21cm cosmic angular power spectrum. The cross-correlation 
can be represented by 

/ dT b 

(T em Tab) = To(l + 6)xn 



dS 



The angular power spectrum weighted by window function W v (r) 

Ccm-ab(^) = To 



and under small-angle approximation is 



^P- (^) [P^Wr)+x R fP M {l/r)]dr. 



(25) 



(26) 



Note that the ionized bubbles surrounding the first galaxies are excised in the above computation to account for the null effect 
of 21cm absorption from neutral hydrogen. 



3.6 Cross-correlation: C m h,_ a b(-Q 

Unlike the above two cross-correlations which may manifest themselves only in the rapid transition phase, the 21cm emission 
of minihalos should naturally exhibit a correlation with the 21cm absorption of uncollapsed IGM through the entire period 
of the Dark Ages till the Grey Ages. The cross-correlation between these two components is 



d 2 n(M, z) 



dMdV 



The fluctuation arises purely from the density perturbation 5, and therefore its power spectrum is simply 
P^)=(^)^±^(^)fPu(k). (28) 



2k B ) H(z) V dS 

Finally, we can estimate the angular power spectrum modulated by window function W v (r) through 

Cmh- 



One may argue that the fraction of neutral hydrogen in the above calculations may be overestimated because the contribution 
of minihalos to the total baryonic matter of the universe has not been excluded in the treatment of the uncollapsed IGM. To 
clarify the point, we have calculated the hydrogen mass inside minihalos at different redshifts and found that its contribution 
to the total fraction of global neutral hydrogen is negligibly small. 



4 RESULTS 

We adopt a frequency resolution of Ais — 0.1MHz to conduct our numerical computations of the angular power spectra of 
three auto-correlations and three cross-correlations. The total 21cm angular power spectrum is then obtained by summing 
up the six components following Eq.(3). We demonstrate all the resultant angular power spectra in the conventional way, 
ST = [£(£+ l)C(£)/27r] 1//2 , which provides a direct measurement of rms fluctuation of the 21cm brightness temperature in the 
corresponding multipole range. 

Figure 1 shows the angular power spectra of three auto-correlations (C om -cm, C m h-mh and C a b-ab ) and three cross- 
correlations (Cem-mh) C om _ab and C m h-ab) calculated at four different redshifts, z=10, 15, 20 and 25. The results of the 
three auto-correlations are certainly not novel and have already been demonstrated in the literature (see, for examples, 
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Figure 1. The predicted angular power spectra of 21cm brightness fluctuations at four different redshifts, (a) z = 10, (b) z = 15, (c) 
z = 20 and (d) z = 25. Three auto-correlations (red lines) and three cross-correlations (green lines) are illustrated by dotted curve 
(Cera-™), dot-dashed curve (C mh - mh ), long dashed curve (C ab - ab ), dotted curve (C em ^ mh ), dot-dashed curve (C mh - ab ), and long 
dashed curve (C em _ a (,), respectively. Solid black curve represents the total angular power spectrum. Note that each of the cross-correlation 
components C(£) has been multiplied by a factor of 2 to account for its effective contribution to the total power spectrum [see Eq.(3)]. 
For the two negative cross-correlation components relevant to the uncollapsed IGM, their absolute values are utilized. 



Iliev et al-lbooi lLoeb fc Zaldarriaeal Eooi : Zaldarriaga et al. 2004). As what is expected, the overall power spectrum for £ 



up to 10 is governed by the auto-correlation term, C C m- C m at lower redshifts and C a t>-ab at higher redshifts, respectively, 
and the transition occurs at 15 < z < 20. The auto-correlation of minihalos contributes only a small fraction to the total 
power spectrum. Nonetheless, the cross-correlation of minihalos with other two mechanisms constitutes the second largest 
components to the total power spectrum: C om _ m h at redshifts below z ~ 15 and C m h- a b at redshifts beyond z ~ 20, though 
their absolute values are still small as compared with the contribution of corresponding auto-correlation (C cm -cm at lower 
redshifts and C a t>-ab at higher redshifts). 

While the effect of cross-correlations on the overall 21cm angular power spectrum is insignificant, the cross-correlation 
may exceed the auto-correlation during the transition phase which occurs at a certain redshift between z w 15 and z w 20 
in terms of the present models. To demonstrate the point clearly, we take the 21cm power spectrum of each component at 
two multipoles, i = 10 3 and £ — 10 4 , and illustrate in Figure 2 and Figure 3 its time evolution as well as its fractions in the 
total power spectra. Note that under certain circumstances some fractions may be larger than unity because of the negative 
contribution of the two terms, C cm - a b and C m h- a b, to the total power, which is actually a good indicator of the significance 
of cross-correlations at the corresponding multipole range. For comparison, we have also displayed the results when only the 
auto-correlations are taken into account. It turns out that inclusion of the cross-correlations leads to a suppression of the 
angular power spectrum during the transition phase, and the extent and position of such a suppression varies with multipole 
£. For the example of i = 10 3 shown in Figure 2, there is a drop of up to an order of magnitude around z = 17 — 18 with 
regard to the angular power spectrum without the inclusion of cross-correlations for the same £. 

It is not difficult to understand the cause for the above power suppression during the transition phase. If we neglect the 
contribution of minihalos, the total power spectrum reduces to C(£) = C cm -cm(£) + C a b- a b(^) — 2|C em _ a b(^)|- Before and 
after the transition, C(£) is dominated by either C em _ cm (^) or C a b- a b(^), depending on redshifts. During the transition phase, 
however, the two effects make more or less the same contribution to the total power spectrum, or the auto-correlation is 
roughly the same as the cross-correlation in amplitude, hence reducing greatly the power. In other words, the two components 
almost cancel each other because the cross-correlation has a negative sign. We illustrate the point in Figure 4 by plotting 
the angular power spectrum at redshift z = 17.5, the transition phase in terms of our present models. It appears that a 
combination of the auto- and cross-correlations yields a much lower angular power spectrum. 
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Figure 2. Redshift dependence of 21cm angular power spectra at multipole I = 10 3 . The left panels (a, c) and right panels (b, d) show 
the results without and with the inclusion of cross-correlations, respectively. The top panels (a, b) are the angular power spectra, while 
the bottom panels (c, d) display the fraction of each component in the total angular power spectra. We use the same notations for all 
the curves as in Fig.l. 
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Figure 3. The same as Fig. 2 but for I = 10 



5 DISCUSSION AND CONCLUSIONS 

Using a simple analytic prescription of both 21cm emission and absorption of neutral hydrogen in the epoch of the Grey Ages, 
we have calculated the cross-correlation between three major mechanisms of cosmic 21 cm radiation, namely, the emission 
from neutral hydrogen surrounding the ionized bubbles of first galaxies, the emission of neutral hydrogen inside minihalos, 
and the absorption of diffuse neutral hydrogen against the CMB radiation. It has been shown that inclusion of the three cross 
correlations does not alter the overall angular power spectrum of cosmic 21cm background. However, during the transition 
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Figure 4. Contribution of each correlation component to the total angular power spectrum at the transition phase of z = 17.5. The 
same notations as in Fie. 2 are used for all the curves. Note again that the absolute values of Ccm— ab 

and C m h- a b are shown. 

of 21cm radiation field from emission-dominated phase to absorption-dominated phase, the cross-correlation may play a key 
role in the determination of the shape and magnitude of the 21cm power spectrum. This arises primarily from the negative 
contribution of cross-correlation between the absorption of diffuse IGM and the other two emission mechanisms, which leads 
to a significant power suppression on the fluctuation of 21cm radiation field during the transition phase. In our naive models 
for the above three mechanisms, especially the oversimplified scenario for the 21cm emission from the IGM associated with 
first galaxies, the transition occurs around redshift z ~ 17 — 18, which varies slightly with multipole I. 

The effect of cross-correlations on the 21cm angular power spectrum could be detected as a remarkable drop of power at 
a certain frequency range corresponding to the transition phase. The precise value of the transition frequency depends on the 
combined effects of various mechanisms of 21cm radiation in the epoch of the Grey Ages. An extension of the present work 
to the adoption of more sophisticated models for the 21cm radiation field at the Grey Ages may provide a more accurate 
prediction of the 21cm angular power spectrum. 

Finally, we must caution that the effect of the cross-correlations on the total cosmic 21cm power spectrum could be 
considerably exaggerated as a result of our 'toy model', and some of our numerical results may even be artificial because of 
the 'artificial' split between the three major mechanisms of 21cm radiation. Numerical simulation should be applied to clarify 
the issue. 
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